Molecular spin-liquid state in spin-3/2 frustrated spinel HgCr 2 04 
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A hexamer-type spin excitation seen in spinel chromates ACr2 04 (A = Mg, Zn, or Cd) is the 
representative spin-liquid-like state caused by geometrical frustration. To clarify an origin of the 
state, we comparatively studied spin excitations in an isomorphic material HgCr204 by inelastic 
neutron scattering, and observed a different molecular-type excitation. Numerical analyses per- 
formed using model Hamiltonians suggest that these two types of spin excitations originate from 
a spin-3/2 molecular singlet hidden in a magnetically ordered phase. The difference between the 
molecular types is explained by the difference in the kind of exchange interactions occurring in the 
chromates. 

PACS numbers: 75.30.-m, 75.40.Gb, 75.50.Xx, 75.50.-y, 78.70.Nx 
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Remarkable developments have been made in the fields 
of solid-state physics in magnetic systems and strongly 
correlated electron systems using the concept of geo- 
metrical spin frustration. In frustrated systems, not 
all classical-spin pairs can be arranged antiferromagneti- 
cally on a triangle lattice and a tetrahedral lattice [l|-|3|- 
Therefore, frustration suppresses magnetic ordering and 
promotes spin-liquid-like fluctuations (zero-energy exci- 
tations) in a low-temperature paramagnetic phase, e.g., 
spin molecules, spin ices, and spin vortices [4-8]. Re- 
cently, dynamical spin molecules were found to exist as 
non-dispersive gapped excitation modes within a mag- 
netically ordered phase, where frustration was assumed 
to be relieved by a lattice distortion [9(. 

Formation of a spin molecule is one of the represen- 
tative effects of geometrical frustration in a system 
A spin molecule refers to a spin cluster that is spa- 
tially confined within a geometrical shape such as an 
atomic molecule. Intermolecular correlation among spin 
molecules is negligible in comparison to the intramolec- 
ular correlation; this behavior is similar to the quasipar- 
ticle approximation in Landau's Fermi liquid theory 
For example, a paramagnetic phase in the spinel ACT2O4 
(nonmagnetic A — Mg, Zn, Cd) is known to undergo non- 
frustrated hexamer-type spin fluctuations so as to avoid 
frustration (hexa-I in Fig.^b)) @,ErIll|. The magnetic 
ions Cr 3+ (fog) 3 , spin S = 3/2) form a corner-sharing 
tetrahedral lattice called a pyrochlore lattice; such a lat- 
tice has kagome and triangle planes stacked alternately 
along the [111] direction and is geometrically frustrated. 

HgCr 2 04 is isomorphic to the aforementioned spinel 
chromates. Its Curie- Weiss temperature &cw has been 
estimated to be —32 K, indicating an antiferromagnetic 
spin correlation 12j. It undergoes antiferromagnetic 



&cw owing to frustration [l2j, |l3j. A complex mag- 
netic structure of HgCr 2 04 determined by powder neu- 
tron diffraction also suggests that antiferromagnetic first- 
neighbor exchange interactions govern the spin system, 
as is the case with the other chromates fl3j [. 

However, the magnitude of &cw and the degree of 
frustration \@cw /Tn\ — 5 are both lower than those of 
the other chromates; for example, the &cw values of Mg, 
Zn, and Cd chromates are about —370, —390, and —70 K, 
and their correspondin g; \& nw /Tn\ values are about 30, 
31, and 9, respectively 12|, LLJ]. A low Qcw indicates a 
scaling-down of exchange interactions. Alow \@cw/Tn\, 
in turn, indicates a suppression of frustration by rela- 
tively strong further-neighbor exchange interactions and 
spin-lattice coupling, which also explain the wide mag- 
netization plateau observed under a high magnetic field 
from 10 to 27 T 0, [I3-[l7 1 . 



long-range ordering with propagation vectors (1/2,0,1) 
and (1,0,0) below Tat = 6 K, which is much lower than 



Thus, HgCr204 might exhibit spin excitations other 
than hexa-I. To the best of our knowledge, no inelas- 
tic neutron scattering experiments on HgCr204 have yet 
been performed, probably because of the strong neutron 
absorption of Hg nuclei. 

In this study, therefore, we performed powder inelas- 
tic neutron scattering experiments on HgCi'204 above 
and below Tjy without a magnetic field. This Letter re- 
ports the discovery of another type of dynamical spin 
molecules from the experiments. It also discusses the 
origin of the difference between molecular types and pro- 
poses a quantum-mechanical picture for spin molecules 
by numerical analyses using model Hamiltonians. 

A powder sample of HgCr204 was synthesized by ther- 
mal decomposition of Hg2Cr04 in an evacuated silica 
tube Some preliminary neutron experiments were 

performed on the triple-axis spectrometer TAS-2 of the 
Japan Atomic Energy Agency (JAEA), installed in the 
thermal guide tube of the JRR-3 reactor in JAEA. In- 
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elastic neutron scattering experiments were performed 
on the triple-axis spectrometer HER (Cl-1) at the Insti- 
tute for Solid State Physics (ISSP), University of Tokyo; 
this spectrometer is installed in the cold guide tube of 
the same reactor. The energy of the final neutrons, Ef, 
was fixed to 3.6 meV with a horizontal collimation se- 
quence of guide-open-radial-open, where the radial colli- 
mator has 3 blank channels. A horizontal focusing ana- 
lyzer merged the scattered neutrons in a range of about 
7 degrees scattering of angle. A cooled Be filter and a 
pyrolytic graphite Bragg-reflection filter efficiently elim- 
inated the half-lambda contamination. The powder (4.5 
g) was filled in a thin aluminum foil, which was then 
shaped into a hollow cylinder of diameter 25 mm in or- 
der to minimize the strong neutron absorption. Then, 
the sample was enclosed in an aluminum container with 
4 He exchange gas, which was placed under the cold head 
of a closed-cycle 4 He refrigerator. 

FiguresQJa) andQJb) show the scattering intensity dis- 
tributions in momentum Q and energy E space, mea- 
sured above and below TV, respectively. Above TV, 
quasielastic scattering with two peaks at Q ~ 0.9 and 1.5 
A -1 is observed, as shown in Fig.[lja). In contrast, below 
Tjy, this double-Q scattering becomes gapped resonance- 
like excitations around 1.5 meV, as shown in Fig. [TJb). 
The central values of resonance energies at the two Q's 
seem to be slightly different. 

The Q dependence of -/^-integrated intensity with 
higher statistics is shown in the top panel of Fig. HJc). A 
homothetic double- Q curve is observed at the two tem- 
peratures. The peak around 1.5 A -1 is attributable to 
hexa-I, since the Q value coincides with that observed for 
MgCr 2 4 and ZnCr 2 4 [1 El]. The emergence of the 
peak at around 0.9 A -1 is characteristic of HgCr204. 

In order to confirm that the resonance-like excitations 
consist of two modes with different Q peak positions and 
similar excitation energies, we compared the Q depen- 
dences of intensities integrated in two narrow E ranges 
(see the middle and bottom panels in Fig. Hfc)). The 
1.7 meV data show only the single peak around 0.9 A -1 
(middle), whereas the 1.1 meV data show the double- 
ts structure with relatively stronger intensity around 1.5 
A -1 (bottom), evidencing the coexistence of the two 
modes. 

Figure [Tf d) shows E spectra integrated from Q = 0.5 
to 2.0 A -1 and measured at several temperatures around 
Tjv- With decreasing temperature, the shape of the 
quasielastic spectra observed above TV changes abruptly 
to the gapped shape below TV, similar to a first-order 
transition. The spectral profile appreciably becomes 
asymmetric at 3 K. A double-Gaussian fitting with a con- 
stant background and an elastic line gives central values 
of 1.3 and 1.7 meV for the two modes. 

To summarize, we discovered two almost degenerated 
spin-excitation modes below TV: the first mode is cen- 
tered around 1.3 meV and 1.5 A -1 and is most probably 
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FIG. 1: (Color online) Measured inelastic neutron scattering 
data for powder HgCr2 04. (a)(b) Scattering intensity dis- 
tributions in (Q,E) space above and below Tjv, respectively. 
The vertical tone indicates the scattering intensity in arbi- 
trary units, (c) Q dependence of ^-integrated intensity. The 
integration ranges are given in each panel. The solid lines are 
a guide to the eye. (d) E dependence of intensity integrated 
from Q = 0.5 to 2.0 A -1 measured at several temperatures. 
In the data for 10, 7, and 5 K, the vertical origins are shifted 
to the height shown by horizontal solid lines on the right side, 
and the solid curves are a guide to the eye. The solid curve 
for the 3 K data is obtained by a multi- Gaussian fitting, indi- 
cated by the dotted and broken lines. The errors are smaller 
than the symbols used. 



based on hexa-I, and the second mode is centered around 
1.7 meV and 0.9 A -1 . The modes are quasielastic above 
TV- 

Next, in order to determine the dynamical spin struc- 
tures, we analyzed the Q dependence of quasielastic exci- 
tations, though it is difficult without single-crystal data. 
The cross section of spin molecules is described by 

5(Q) = C 1 |F(Q)| 2 ^^exp( i Q.r J )| 2 ^ +B.G, (1) 

where C\ is a scale factor of the experimental scatter- 
ing intensity, B.G. is a constant background, F(Q) is 
the magnetic form factor of Cr 3+ , for which the Watson- 
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Freeman form factor was used below [19j . j denotes the 
site of Cr 3+ , N is the total number of sites in a molecule, 
Tj is the position, and Sj takes only ±1 corresponding to 
collinear spins that dynamically fluctuate in arbitrary di- 
rections [a, Q . The parentheses indicate an orientational 
average with resolution convolution. Using Eq. ([TJ, we 
determined a kagome-star dodecamer model, which can 
be regarded as a superposition of the first-neighbor hex- 
amer (hexa-I) and a second- neighbor hexamer (hexa-II), 
as shown in Fig. EJb) . The molecules are formed on the 
kagome planes. Figure [^a) shows the calculated Q de- 
pendence of intensity for the dodecamer (solid curve), 
which is in good agreement with the experimental data. 

We also note that the dodecamer model cannot be 
definitively distinguished from the model of the superpo- 
sition of individual hexa-I and hexa-II. Figure [Ha) also 
shows the calculated curves for hexa-I (broken) and hexa- 
II (dotted) , whose summation with equal weight gives the 
same curve as the solid one. Hereafter, we refer to both 
these models as the dodecamer model without any dis- 
tinction. 

Using the dodecamer model, we analyzed the 
resonance-like excitations observed below TV by assum- 
ing the following effective Hamiltonians: 



-ff hexa-I — J\ y Si ■ Sj, 
-^hexa-II = J2 J]] Si ' Sj, 



(2) 
(3) 



where 5, = 3/2, i and j denote the hexagonal sites 1 
to 6, J2{ij)-i an d Y^,(ij)-n denote summation over all 
first-neighbor and second-neighbor S pairs (not dou- 
bly counted), and J\ and J2 denote the first- and 
second-neighbor exchange interactions that are antifer- 
romagnetic, as expected from the quasielastic excita- 
tions (Fig. Mh)). We used 4096 (= 4 6 ) basis states of 
\S{ , SI, SI, SI SI, SI), where Sf = ±3/2 and ±1/2. Af- 
ter obtaining molecular ground states | Ao) and nth ex- 
cited states I A n ) with excitation energy E n by exact 
diagonalization of Eqs. ((2} and ([3|), the cross section of 
magnetic inelastic neutron scattering at E — E n can be 
calculated by 



S(Q,E n ) = C 2 \F(Q)\ 2 5(hw-E n )( £ (<^ 



va,/3=l 



QaQ/3 . 

' IQI 2 ' 



JV \ 

(Ao I Sf I A„)(A„ I §f, I A ) exp{zQ • (r, - ry)} 
ij'=i / 

+ B.G, (4) 

where C2 is a scale factor of the inelastic scattering in- 
tensity, a and /? are the components of the Cartesian 
coordinate, and S is a spin operator 2l|, |22j. The di- 
rectional term (8 a p — Q a Qp/\Q\ 2 ) can be considered as 
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FIG. 2: (Color online) Dynamical structure modeling of spin- 
excitation modes in HgC^CU. (a) Experimental Q depen- 
dence (symbols) and calculated Q dependence (curves). The 
experimental data are obtained by subtracting a low-Q back- 
ground component fitted with a Lorentzian from the raw 
data (top panel in Fig. 0Jc)). The latter curves are obtained 
from a kagome-star dodecamer model consisting of hexa-I and 
hexa-II (b) , and show the common results calculated indepen- 
dently using Eqs. JTJ and Q (see the text). The broken and 
dotted curves broaden approximately with reciprocals of the 
hexa-I and hexa-II spatial sizes around the main peaks at 
Q ~ 1.5 and 0.9 A -1 , respectively, (b) Kagome-star dode- 
camer model. The subunits hexa-I and hexa-II are depicted 
by bold lines. The closed and open circles indicate up and 
down spins, respectively, fluctuating dynamically in arbitrary 
directions. Exchange interactions Ji and J2 are also defined, 
(c) Energy schemes of hexa-I and hexa-II obtained by ex- 
act diagonalization of Eqs. and ©. The hexa-I part in 
HgCr2C>4 is converted into hexa-I in the other chromates by 
tuning the energy scale [2p( ] . 



a constant by taking the orientational average over dy- 
namically fluctuating molecules. 

Figure Etc) shows energy schemes obtained from 
Eqs. and §3§ with antiferromagnetic J\ = 1.8 meV 
and J2 =2.4 meV; these are in good agreement with the 
experimental resonance energies. Interestingly, the two 
spin-3/2 hexamer systems have a nonmagnetic singlet 
ground state (S tot — 0) and triplet first-excited states 
(S tot = 1), where S tot = £? =1 Si. The states are repre- 
sented by complex linear combinations of S z = ±3/2 and 
±1/2. A direct product of the hexa-I and hexa-II singlets 
generates a dodecamer-singlet ground state, and the two 
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types of hexamer singlet-triplet excitations correspond 
to its intra-activations. A similar dodecamer singlet unit 
with six spin-1/2 dimers, called a pinwheel, was recently 
discovered in a kagome material by single-crystal inelastic 
neutron scattering [23J. The calculated Q dependences 
of intensity for the two types of hexamer excitations, ob- 
tained from Eq. (|4|, are completely nomothetic to those 
for the quasielastic excitations (broken and dotted curves 
in Fig. Ufa)), respectively, whose summation is again in 
good agreement with the experimental data. 

For the hexa-I part, we also calculated the inten- 
sity distributions in the hkO and hhl zones from both 
Eqs. (HJ) and §4$ without powder orientational averag- 
ing. The patterns obtained from both these equations 
are the same, which is consistent with the single- crystal 
data measured below and above Tm in MgCr204 0- 

In this way, spin excitations above and below TV in 
HgCr2C>4 were explained by the kagome-star dodecamer 
model with hexa-I and hexa-II, which is based on the 
proximity of J\ and J2. Then, the following question 
arises: what is the reason for this proximity in HgCr204? 
In chromates, J\ involves a direct exchange with a Cr-Cr 
direct overlap of t^ g orbitals, and Ji involves superex- 
change with a Cr-O-A-O-Cr path, with a Cr-O-A bond 
angle of about 125°. The magnitude of J\ weakens with 
increasing Cr-Cr distance, and that of J2 strengthens 
with increasing bond angles up to 180° according to the 
Goodenough-Kanamori rules. Superexchange also occurs 
with 90° Cr-O-Cr paths, which can be neglected accord- 
ing to the rules. Meanwhile, as the ionic radius of an A 2+ 
cation increases from Mg 2+ to Hg 2+ , the Cr-Cr distance 
increases from 2.94 A to 3.06 A and the Cr-O-A angle in- 
creases from 122° to 131° [H,[25j]. Thus, Ji is expected 
to weaken and J2 is expected to strengthen in HgCr204. 
This tendency is also supported by recent LSDA+C7 cal- 
culations (2^ . 

The above analyses revealed that the singlet-triplet 
picture, which is usually based on the quantum spin-1/2 
dimer, can be extended to spin-3/2 hexamers and the 
dodecamer. Singlet formation does not imply complete 
disappearance of magnetic moments and can coexist with 
magnetic order, since the g factor is arbitrary in the anal- 
yses. In fact, powder neutron diffraction experiments 
revealed that ^4Cr204 (A = Mg, Zn, Hg) exhibits an or- 
dered moment with only 2.2, 2.0, and 1.7 /is 13, 27, 28j |. 
These values are around 1 /ib lower than the full value of 
3 /ie, which is consistent with the partial singlet forma- 
tion conjecture. Thus, in the sense of the singlet, mag- 
netically ordered phases in ^4Cr204 could be cited as a 
new class of quantum spin liquid with classical magnetic 
long-range order. Then, quasielastic spin excitations in 
a paramagnetic phase, whose the origin has been one of 
the longstanding issues related to frustration, could be 
interpreted as its thermally fluctuating spin-liquid state. 

In summary, we discovered almost degenerated spin- 
excitation modes in HgCr204 by powder inelastic neu- 



tron scattering. The Q dependence of intensity demon- 
strates the coexistence of the normal hexa-I-type exci- 
tation and another type of excitation. On the basis of 
dynamical spin structure analyses, we propose the use of 
the kagome-star dodecamer model with hexa-I and hexa- 
II, which is based on the proximity between J\ and J2. 
Further, by numerically analyzing model Hamiltonians, 
we proposed a quantum-mechanical picture that is ap- 
plicable to all chromates: the quantum spin liquid (spin- 
3/2 hexamer/dodecamer singlet) coexists with magnetic 
long-range order as ground states below Tjy, and the 
quasielastic excitations correspond to its thermally fluc- 
tuating spin- liquid state above Tn ■ Further experimental 
and theoretical studies will be needed to prove the spin 
liquid conjecture. 
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